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We present a direct calculation of the cross section for 
the reaction 3 He(e, e'p) including the radiation tail orig- 
inating from bremsstrahlung processes. This calculation 
is compared to measured cross sections. The calculation 
is carried out from within a Monte Carlo simulation pro- 
gram so that acceptance-averaging effects, along with a 
subset of possible energy losses, are taken into account. 
Excellent agreement is obtained between our calculation 
and measured data, after a correction factor for higher- 
order bremsstrahlung is devised and applied to the tail. 
Industry-standard radiative corrections fail miserably for 
these data, and we use the results of our calculation to 
dissect the failure. Implications for design and analysis 
of experiments in the Jefferson-Lab energy domain are 
discussed. 
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I. INTRODUCTION 

The application of what are commonly called radiative 
corrections is an integral part of doing nuclear physics 
with beams of electrons. In an electron-scattering ex- 
periment, the probe is considered to be a virtual photon. 
This photon is exchanged between a beam electron and a 
target nucleus, thereby transferring energy and momen- 
tum to the target from the the electron, which is thereby 
scattered. Unfortunately, these electrons also copiously 
emit real photons which are not normally observed in ex- 
periments. Thus, either the theoretical calculations with 
which data are compared must include these processes 
(and they normally do not), or the data must somehow 
be corrected for these effects so that they can be com- 
pared to calculations which are based on single virtual- 
photon exchange. The standard choice is to "radiatively 
unfold" the experimental data, which generates a "cor- 
rected spectrum" that can be compared to theoretical 
calculations. 

This paper reports on a study of how these real-photon 
processes affect measurements of (e, e'p) reactions on 



atomic nuclei. Our calculation takes the second ap- 
proach, which is to radiatively correct a theoretical calcu- 
lation so that it can be directly compared to uncorrected 
data. We compare a direct computation of a cross sec- 
tion, including the effects of photon emission, to a specific 
measurement of an (e, e'p) cross section [|l|]2J]. The results 
of applying the standard "radiative unfolding" procedure 
mentioned above to these data are also presented and dis- 
cussed. The comparisons to this particular data set are 
unique in two ways: 

1. the data appear to be well described in the plane- 
wave impulse approximation (PWIA), and 

2. over most of the kinematic range of the measure- 
ment, the real-photon "radiative processes" are re- 
sponsible for nearly the entire cross section. 

The excellent PWIA reproduction of the cross section 
chosen for this study allowed us to use the PWIA in 
carrying out the complex calculations including radiative 
effects, enormously simplifying the task. The dominance 
of radiative strength enables us to make a true test of 
the real-photon emission model without worrying about 
accurately removing physical backgrounds. 

This study is timely for several reasons. First, ex- 
isting procedures for radiative corrections to data have 
been developed for experiments at relatively low (< 500 
MeV) electron beam energy. Refinements or overhauls 
of the procedure may be necessary to apply corrections 
for (e, e'p) experiments with higher-energy beams. The 
experiment studied here was carried out with a beam en- 
ergy of 855 MeV, which bridges the gap between the en- 
ergy domain studied by the labs active in the last decade 
(0.2-0.9 GeV) and the Jefferson Lab energy domain (0.8- 
6.0 GeV). Furthermore, a new class of experiments at 
Jefferson Lab begins to study (e, e'p) reactions in a kine- 
matic domain where the cross sections are expected to be 
small and broadly peaked; radiative strength can easily 
swamp the "true" cross section in these cases. The design 
and analysis of these experiments should make careful 
studies of the radiative contributions to measured cross 



sections. Indeed, such an analysis was the genesis for the 
current work. 

Finally, it became clear to us during the course of 
the project described here that the standard radiative- 
unfolding procedure used for the last decade is of an ad 
hoc nature; it is not based on rigorous theoretical ar- 
guments. We could only find one article published in 
a refereed journal |3(] which specifically addressed radia- 
tive corrections for (e, e'p) reactions. This publication 
is either unknown to most experimentalists, or has been 
ignored for some reason, as a literature search uncov- 
ered only one reference Q to this article, which argued 
that the corrections proposed in p] were impractical since 
they make different assumptions about hadronic portions 
of the corrections, rendering data so corrected inconsis- 
tent with the world proton form-factor data. 

All remaining works we could find addressing radiative 
corrections to experimental data were Ph.D. theses. Es- 
sentially all these works quote the Ph.D. thesis of Quint 
P| as the primary reference. This thesis in turn quotes 
lecture notes of Penner [pi from a summer-school pro- 
ceedings as a primary source, where radiative corrections 
for inclusive (e, e') reactions are discussed. These notes 
clearly state that the correction should be viewed as ap- 
proximate; for example, they recommend an empirical 
adjustment of the calculated tail to give the best fit to 
the data, in cases where the radiation tail dominates the 
cross section. Aside from this problem and possible prob- 
lems in adapting a formalism for (e, e') to correct (e, e'p) 
experiments, we have uncovered several questionable as- 
sumptions in this standard procedure, which we address 
in this article. 

By contrast, we base our calculations on a published 

first-order QED calculation for the radiation-tail cross 
section. Our work extends their result to kinematically- 
complete reactions and to higher-order bremsstrahlung 
radiation. 

We are unaware of any previously-published similar 
study. We hope this paper will give some indication of 
how urgently new theoretical work is needed, and in what 
directions that work should proceed. Since the topic of 
radiative corrections is often viewed as an arcane subject 
which is best avoided, we present in the following sec- 
tions a review of the relevant electromagnetic processes, 
a review of (e, e'p) phenomenology, and an explanation 
on how radiative effects distort (e, e'p) reaction data. 



II. REVIEW OF RADIATIVE PROCESSES IN 
ELECTRON SCATTERING 

This section presents a review of the most important 
processes via which electrons emit real photons during 
interactions with nuclei. We emphasize here that this 
section is a review, meant to place these processes in the 
context of the reaction we study. Much of the conceptual 
work here, and formal work on bremsstrahlung presented 
in sec. [V], can be found in classic articles p-|lO|]. 



I 




FIG. 1. Schematic diagram of the (e, e'p) reaction (in this 
case off an A = 4 nucleus) to lowest electromagnetic order. An 
incident electron scatters from a target nucleus by exchange 
of a virtual photon, and a target proton is knocked out in the 
process. Symbols next to the various lines show the names 
given the four- momenta for each particle. 




FIG. 2. Example Feynman diagram of bremsstrahlung in 
e-p scattering. A real photon is emitted from the outgoing 
electron. There are three other such diagrams, one each for 
the two proton legs and one for the incident electron leg. 

Fig. [y is a schematic diagram of the (e, e'p) process to 
leading order in the electromagnetic coupling constant a. 
It corresponds to the mental picture usually employed 
by an experimentalist designing or analyzing an experi- 
ment, since it probes the "signal" the experimenter usu- 
ally wants to measure. It also corresponds to the usual 
plane-wave impulse approximation (PWIA) for (e, e'p) 
reactions. For the purpose of the study presented here, 
we have chosen a measurement on 3 He which was per- 
formed in kinematics specifically chosen to optimize the 
accuracy of the PWIA. However, even in the limit that 
the PWIA holds for the hadronic portion of the process, 
the neglect of real-photon emission limits the accuracy 
of PWIA cross sections to at best 20% for electron ener- 
gies above a few hundred MeV. An example diagram for 
this bremsstrahlung process is shown in Fig. |2j. Here the 
spectator nucleons have been omitted from the picture 
for simplicity. 

The process in Fig. is called internal bremsstrahlung 
since it occurs during the (e, e'p) reaction. A similar 
process (termed external bremsstrahlung) takes place in 
the Coulomb fields of other atoms in the target. A related 
process is not particularly relevant for understanding the 
effect of radiative processes on the (e, e'p) spectra, but 




FIG. 3. Example diagram for the virtual photon correction. 

must be included in any consistent calculation. This is 
the process in which two virtual photons are emitted, and 
an example diagram is shown in Fig. 0. Such diagrams 
are generally termed "virtual photon corrections." 

For fixed values of the four- momenta e and e', we 
see that the value of the four-momentum transfer q is 
changed in the diagram of Fig. [2| with respect to the 
leading-order process in Fig. [l]. This in general leads to 
a change in the magnitude of the associated cross section 
(as does the vertex renormalization in Fig. m. The extra 
emitted particle in Fig. g leads to a change in the asymp- 
totic kinematics of the reaction as well. This creates an 
ambiguity; for a given measured event, it is impossible to 
tell whether the observed kinematics correspond to those 
of the reaction vertex, or to a different reaction-vertex 
situation accompanied by real photon emission. 

We begin our discussion of this problem by summariz- 
ing the phenomenology of the (e, e'p) reaction in the next 
section. The section afterwards will give a quantitative 
description of the kinematic distortion due to the contri- 
butions of the diagram in Fig. g. We present a rather 
complete summary of the kinematics, since discrepant 
conventions exist in the literature. In our discussions 
below, unless otherwise stated we use the following kine- 
matic conventions: 

• the four-vectors are denoted by the standard sym- 
bol for the associated particle, e.g. A for the target 
four-momentum and p' for that of the knocked-out 
proton. 

• Eg and p g refer to the relativistic energy and three- 
momentum of particle g, thus the four momentum 
for g is, g — (E g ,p g ). The magnitude of the three 
momentum is p g = |p s |. 



III. REVIEW OF (e,e'p) PHENOMENOLOGY 

One of the main reasons why the (e, e'p) reaction has 
been so useful in nuclear physics is that it probes the 
properties of individual nuclear protons in a fairly direct 
manner. The measured particle momenta can be used 
to determine the energy and momentum that the struck 



proton had before the interaction. The only major as- 
sumption involved is that the interactions of the recoil- 
ing (^4— 1) system and knocked-out proton are neglected 
(PWIA) . While the PWIA is not sufficiently accurate for 
a quantitative analysis of (e, e'p) experiments, many ex- 
periments have shown |ll]] that the essential features are 
reasonably preserved and that a straightforward analysis 
is possible. 

The probability of finding the "struck" proton of Fig. [j] 
(to which the four-momentum q is transferred) is a func- 
tion of two parameters: an energy (which we refer to as 
e) and a momentum (to which we refer as p). Various 
equivalent conventions for e exist; we use it to refer to the 
energy necessary to remove the proton from the nucleus. 
p refers to the momentum of the proton relative to the 
nuclear rest frame. 

In our convention, e consists of two parts: e = S p + 
E^ . S p is the proton separation energy for the nucleus 
being bombarded, and E^ is the excitation energy of the 
residual system "i?" of A — 1 nucleons. 

Assuming that the PWIA holds, e and p can be com- 
puted from the kinematic variables measured in (e, e'p) 
experiments. We begin by constructing a four-vector re- 
lation for the process (using the notation of Fig. fil): 



A = 



R + p'. 



(1) 



Assuming the reaction is carried out with a known beam 
energy, a fixed, pure target, and that the four- momenta of 
the scattered electron and knocked-out proton are mea- 
sured in detectors, the kinematics are uniquely deter- 
mined. 



R = (E R ,p R ) = (e-e') + A-p'. 



(2) 



The invariant mass of the (A — 1) system, vi? 2 , yields 
m R ; an experimental missing energy?] is computed as 



E m = m R + m p - tua- 



(3) 



When PWIA holds, e = E m . Similarly, an experimental 
missing momentum is defined as 



PR 



1 Some authors use E m to denote the unmeasured "missing" 
energy in the reaction, which thus includes the kinetic energy 
of the recoiling undetected system. This terminology is his- 
torically correct, since in early experiments with low-energy 
beams on heavy targets, the recoil kinetic energy was negligi- 
ble. E m became synonymous with the binding energy. Later, 
approximate corrections were used to remove the recoil en- 
ergy. The use of relativistic invariants eliminates the need for 
approximations. Our value might be more properly termed 
the "missing mass" since all other forms of energy have been 
accounted for; nevertheless, we stick with the historical term. 



When PWIA holds, the residual system is a spectator 
and thus must have had the same momentum p m before 
the interaction. Since the nucleus as a whole was initially 
at rest, p = -p m = -p R . 



IV. THE KINEMATICAL EFFECTS OF 
BREMSSTRAHLUNG 

If we add a real photon as shown in Fig. to one of 
the external legs in Fig. |l|, we must account for it in the 
four-momentum conservation relation. Here we keep us- 
ing E m andp m for the names of the measured quantities. 
We indicate by use of the extra subscript v (e.g., p m ,v) 
the corresponding quantity at the qA (virtual-photon- 
nucleus) reaction vertex in the case that the actual reac- 
tion involved emission of a real photon. 

The four-momentum conservation relation becomes 



R = (e 



A - p' - 7, 



(4) 



where 7 = (fc,k) refers to the real photon's four- 
momentum. The three-vector component of this equa- 
tion yields 



Pb.,v = q — P P ' — k, so that 

Pm,t; Pm ^ • 



(5) 
(6) 



Thus the deduced missing momentum is p m = p m ,v + 
k. The effect of the real photon emission on E m is not 
obvious when using four-momentum algebra; instead, we 
note that the zeroth component of fl) leads to 

m R , v + ni p -m A = (E e - E e > - T p > - T R>V ) - k, (7) 

where T refers to the particle's kinetic energy and k is the 
photon energy. The left-hand side of this relation is E m ^ v , 
and the value in parentheses is, to a good approximation, 
what one would deduce for E m if one is ignorant of the 
photon emission. Thus E m ss £ m> „ + k. The relation is 
approximate since one measures T R , not T RiV . However, 
the difference is in most cases quite small. 

The bremsstrahlung-photon emission thus causes 
cross-section strength which would normally populate 
(E m ,v,Pm.v) to instead be redistributed over a range of 
values (E m ,p m ). Furthermore, the magnitude of this re- 
distributed cross section will be modified since the mo- 
mentum transfer q will be changed. The redistributed 
strength is the origin of the well-known "radiation tail" 
of electron-scattering experiments. In general, the miss- 
ing energy is simply increased by an amount equal to 
the radiated-photon energy. The missing momentum is 
shifted in a kinematic-dependent manner; the relative 
orientation of its vertex value p m ,v and that of the radi- 
ated photon k plays an important role. Fig. [| illustrates 
how measured strength in a particular region of (E m ,p m ) 
is fed, through the bremsstrahlung process, by various re- 
gions of (£,„,„,?„,,„). 



A procedure (see Ref. Q for a good discussion) has 
been developed for "radiatively correcting" (e, e'p) spec- 
tra. The procedure relies on the fact that in reactions for 
which E m is a minimum, there is no tail, only a reduction 
in the cross section due to the absence of that strength 
which has been moved into the tail at larger values of E m . 
The procedure corrects cross section data by beginning 
with this minimum- E m bin, using the Schwinger correc- 
tion to correct for the amount which has been lost to 
the tail. Then the tail itself is computed from this effec- 
tively "unradiated" cross section. The tail contribution 
from this bin is subtracted from all bins at larger miss- 
ing energies. This procedure can be iterated bin-by-bin, 
moving from small E m to large, to remove the radiation 
tail. 

Estimated uncertainties in the value of the correction 
are usually around 10-20%, which is acceptable when the 
correction itself is small. However, for experiments inves- 
tigating the large- E m continuum cross section, the cor- 
rection can become rather large, or the radiation tail can 
even dominate the cross section. In these cases, even a 
10% uncertainty in the corrections can lead to essentially 
zero knowledge of the true "unradiated" cross section. It 
is important to be able to reliably estimate the strength 
of the radiation tail, so cases like this can be avoided 
during the planning stage of an experiment. In the fol- 
lowing section, we describe the procedure for computing 
the cross section, including the radiation tail. 



V. COMPUTATION OF (e, e'p) CROSS SECTIONS 
INCLUDING RADIATIVE PROCESSES 

Our approach to direct computation of the (e, e'p) 
cross-section spectra, including bremsstrahlung pro- 
cesses, is based on the PWIA for the (e, e'p) reaction. 
The PWIA assumption is not a necessary one; however, 
a computation involving a more complete theory would 
be computationally much more intensive. The use of 
PWIA is well-motivated in our case since it works well 
for 3 He(e,e'p), apart from an overall scaling factor; this 
has been observed in other experiments as well (see e.g. 
fll3|). The basic program is as follows: a Monte-Carlo 
simulation is performed in which the four-momenta of 
the scattered electron and knocked-out proton are sam- 
pled over their respective acceptances. The kinemat- 
ics for the (e, e'p) reaction vertex are then modified for 
bremsstrahlung processes according to the corresponding 
probability distributions. PWIA is finally used to com- 
pute the cross section, and any relevant jacobian factors 
are applied. 

The data generated can then be used to form cross- 
section spectra using the same procedures applied when 
analyzing the experimental data. The use of Monte-Carlo 
simulation is important for an additional reason: the 
cross sections are often rapidly varying over the experi- 
mental acceptances. A calculation of cross-section spec- 
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FIG. 4. Example bremsstrahlung trajectories for 3 He(e,e'p) in the (E m ,v,Pm,v) plane for the kinematic settings given in 
Table |. Each line represents a trajectory through which strength from lower values of E m>v can feed into the selected bin of 
measured (E m ,p m ). The four lines show paths for radiation occurring either on the incoming or outgoing electron leg, and also 
for two values of the deduced measured angle between p m and q, Pm - q . 



tra using only the kinematics corresponding to the cen- 
ters of all the acceptances does not usually reproduce the 
experimental spectra. Indeed, when radiative corrections 
are applied to "deradiate" experimental data, often one 
of the biggest uncertainties stems from the acceptance- 
averaging assumption made. We will discuss this point 
in more detail later in the article. Our Monte-Carlo pro- 
cedure allows for acceptance-averaging identical to that 
of the experiment. 



A. PWIA Cross Sections 

The cross section for continuum (e, e'p) reactions, in 
the case where the residual (A — 1) system has a contin- 
uum of possible invariant masses ttlr is given by 



d 6 a 



dfl e i duodVLpi dE p i 



PpiEpi<j e pb>yE m ,p m ). 



(8) 



u> is the electron energy loss, given in the laboratory sys- 
tem by us = E e — E e i ; equivalently it is the zeroth com- 
ponent of the four- vector momentum transfer q — (w, q). 
a ep is the elementary cross section for scattering of an 
electron from a moving nucleon. We used the "ecl" pre- 
scription pj] for this cross section. The differences be- 
tween "ecl" and "cc2" at these kinematics is everywhere 
less than 3 % 111 . Our calculation uses the Simon param- 
eterization ]15|f of the nucleon electromagnetic form fac- 
tors in the computation of a ep - S(E m ,p m ) is the proton 
spectral function, which gives the probability of finding 
a proton in the nucleus with momentum p m and removal 
energy E m . 

The simulation samples particle momenta instead of 
energies, so the cross section must be adjusted before 
direct use by the program. For the continuum case, we 
have to apply the transformations u> =>■ p e > and E p i =>■ 



p p i . The resulting cross section used in the simulation is 
related to that of Eq. (@) by 



d 6 c 



P P > 



d 6 cr 



dVL e idp e i dQ.pi dp p > E 1 dVL e i dujdQ.pi dE p 



(9) 



In the case of reactions leaving the residual system in 
a discrete state, the cross section is given by 



d 5 a 



dQ e idQptdp e 



Pp'Ep'Vep n(pm) 

R 



(10) 



For these reactions E m has a definite value E a , so 
S(E m ,p m ) is replaced by the momentum distribution 
n{p m ), where 



S(E m ,p m ) = n(p m )S(E n 



E n 



(11) 



R is the "recoil factor" (really a jacobian factor trans- 
forming E^ =>■ E m ), and is given by 



R=l 



Ep' Pr ■ Pp> 

Er |p„,| 2 



(12) 



Here we have not included the extra subscript v on 
the kinematic quantities, but it should be understood 
when evaluating the radiative cross sections later in this 
section, that the hadronic cross section terms must be 
evaluated at the hadronic- vertex kinematics. Subscripts 
have been added in that section as a reminder. 

So far, we have only studied light nuclei with only one 
possible discrete transition (the A — 1 ground state), so 
our complete simulation consists of a sum of one discrete 
simulation and one continuum simulation. These two 
simulations are each themselves composed of two sim- 
ulations to handle the different pieces of the radiation 
tail. The procedure is straightforward to extend to more 
complicated situations including additional bound-state 
channels. 



B. External Bremsstrahlung 

External bremsstrahlung is relatively simple to in- 
clude. The simulation program contains a model for the 
actual reaction target. For each event, an interaction 
vertex is chosen randomly within the intersection of the 
beam and target volumes. After this choice, the amount 
of target material traversed by the electron before and 
after the reaction is computed. The electron is trans- 
ported through each material (e.g., target cell walls or 
target gas). After each traversal, a new electron energy 
i s g enerated directly from sampling Tsai's distribution 
|[10f for bremsstrahlung from a thin radiator: 



Iexk( E o,k,t) 



bt 



r(i + bt) 



k 



1 / k 3 / k N 

k I ~^ + 4 \E^, 



(13) 



Here k is the radiated photon energy (or energy lost by 
the electron) , Eq is the energy of the electron upon en- 
tering the radiator material, t is the thickness of the 
radiator material in radiation-length units, b is Tsai's 
bremsstrahlung parameter (see Eq. (4.3) of |10)), and T 
is the usual gamma function. 



C. Internal Bremsstrahlung 

Internal bremsstrahlung is included using the cross sec- 
tions for first-order photon emission derived by Borie and 
Drechscl W. Their derivation made use of the peaking 
approximation, which assumes that bremsstrahlung pho- 
tons are only emitted along either the incident beam di- 
rection, or the direction of the scattered electron's mo- 
mentum. A critical review of the peaking approximation 
can be found in ||. We use their results to estimate the 
validity of the peaking approximation for our kinematics, 
and find that it should be accurate to better than 1%. We 
note here that it is difficult to make blanket statements 
about the peaking approximation, except that it becomes 
increasingly worse for larger bremsstrahlung-photon en- 
ergies. 

The Borie-Drechsel cross section was also derived 
specifically for (e, e'p) reactions to the (A— 1) continuum. 
Part of the present work is an extension of that formalism 
to processes in which the (A — 1) system is in a discrete 
state. We also present a derivation of a correction fac- 
tor which accounts for higher-order bremsstrahlung pro- 
cesses. 

For the continuum case, there is complete kinemat- 
ical freedom for all particles, as long as the invariant 
mass of the (A — 1) system is large enough to be above 
the particle-emission threshold of the (A — 1) nucleus. 
The simulation then samples all kinematic variables (the 
scattered-electron three-momentum, the ejected proton 



three-momentum, and the emitted photon momentum). 
In this case, the relevant cross section is given by Ref. J7), 
but we repeat it here with different notation and in a 
form consistent with our results for the discrete case. 
Unless otherwise specified, all kinematic quantities re- 
fer to the asymptotic situation, i.e., what would be as- 
signed if one was not aware that a real photon had been 
emitted. These cross sections have a two-term structure 
which arises from the peaking approximation. The first 
term below corresponds to "preradiation" (photon emis- 
sion along the beam direction) and the second term cor- 
responds to "postradiation" . 



d r a 



dkdiodpp' dQ e / dVL p i 



Im A\dk) t 

' dr' 



d 6 (e- 7 ,e') 



\dk)e> 



dajydpp'dfl 
d 6 (e,e' + 7 ) 



e' ,v^ ^p' 



duj v dp p ' d£l e > iV dTt 



-}. (14) 



The factor / mp corresponds to our "multi-photon" cor- 
rection factor which we will discuss later; / mp = 1 cor- 
responds to the result published in JTj. In this section, 
we can use k for both the energy and the momentum of 
7 since it is a real photon. The dr/dk terms are essen- 
tially jacobian factors for the photon emission. They will 
be given below. Finally, the two cross sections on the 
right-hand side of Eq. (nj) arc the usual "unradiated" 
cross sections, and must be evaluated at the vertex val- 
ues. This is why for example the first cross section is a 
function of e — 7 (the beam four-momentum adjusted for 
photon emission before the interaction) rather than of e. 
For the discrete case, the kinematics are over- 
determined. Since both the scattered electron and 
ejected proton are "detected" in the simulation, but 
the photon is not, we sample over the six-dimensional 
(Pe',Pp') space. For each point in this space, photon en- 
ergies can be chosen which belong to this coordinate and 
as well result in the correct invariant mass of the (A — 1) 
system. k e is the real-photon energy in the case that 
the photon is emitted along the direction of the incident 
electron, and k e j is that for the case of photon emission 
along the scattered-electron direction. These values are 
in general not the same (as opposed to the continuum 
case of Eq. dlj), where the values of k were the same). 



A 2 



' 2p e ' (jV + Pe' - p e /2) - (gy + p e 
2[A + p e • (p p / +p e < ~Pe)] 



tVpJ — 



A 2 - 



(15) 



- 2p e , • (p p / + p e //2 - p e ) - (p p / - p e f 



2[A + p e ' • (p p > +p e > ~Pe)] 

A = tjia + u> — E p i 

_ Pe.e' 
Pe,e' — I 1 

|Pe,e'| 



(16) 
(17) 

(18) 



Again here u> refers to the observed, asymptotic value, not 
that at the vertex. Similarly, A includes the total energy 
of both the recoiling hadronic system and the radiated 
photon since u> is the asymptotic value. 
The associated cross section is 



d 6 a 



duidppi dQ e i dflpi 



dk (dr\ d 5 c(e - 7e,e') 
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/ mp has the same meaning as in the preceding paragraph. 
7e, e ' are the four-momenta corresponding to k eie i. The 
above result (setting / mp = 1 for the moment) was gen- 
erated by substituting Eq. (g), coupled with the discrete- 
final-state expression for the spectral function (Eq. (pd|)), 
into the Boric-Drechsel formula for the radiation tail 
(Eq. (|14|)). The integral over dE mv was formally car- 
ried out by converting it, with the help of appropriate 
jacobian factors, to an integral over dk. The kinemat- 
ical factors for photon emission, one each for pre- and 
postradiation, are given by 
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In the continuum case, k e and k e i are identical (the sam- 
pled photon energy). The jacobian factors transforming 
the cross section from differential in k e , e ' to differential 
in p p i are 
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D. Schwinger Correction 

The internal-bremsstrahlung cross section given above 
becomes singular as the radiated-photon energy goes to 
zero. Hence it can not be used to provide the complete 
radiated cross section. The classic technique is to choose 
a cutoff energy AE which is comparable to the exper- 
imental energy resolution; a radiation tail is generated 
with photon energies between AE and the full energy of 
the radiating electron. The remaining cross section in the 
originating kinematic bin {i.e., the cross section for this 
particular reaction where the total energy radiated away 
by real photons is less than AE) is calculated by comput- 
ing the cross section without the internal bremsstrahlung 
graphs, and then reducing this cross section to account 
for that strength which was moved into the radiation 
tail. This reduction factor is called the Schwinger cor- 
rection. For brevity, in discussions below we will refer to 
the strength remaining in the original kinematic bin as 
the "unradiated strength" . 

The formalism we use for the Schwinger correction is 
due to Penner |6|, and is written as 
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where S r is the first-order correction for internal 
bremsstrahlung. Penner's formulation is based on that 
of Maximon (the expression at the bottom of p. 199 
of ref. [pi]) with the addition of kinematic recoil cor- 
rections proposed by Tsai |17[ |. Furthermore, the part 
of this correction corresponding to real-photon emission 
{S r ) has been exponentiated. Exponentiation of this first- 
order correction was suggested by Schwinger || as a 
means of accounting for higher-order (multiple-photon) 
bremsstrahlung. S v is the correction for virtual-photon 
loops at the reaction vertex. The two 5 factors are given 
by 
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with the "recoil factors" given by 

C, — 1 H r~0- ~ cos# e ') and 
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Li2{x) is the Spence function defined by 

ln(l - y) 



L 2 {x) = - 



V 



-dy. 



Finally, Q 2 denotes the standard square of the four- 
momentum transfer, Q 2 — — (e — e') . 



This version of the Schwinger correction does not ac- 
count for possible real-photon emission by the hadrons in- 
volved in the reaction. Makins pj has noted that this pro- 
cess may begin to become important momentum trans- 
fers Q 2 > 1 GcV/c. Penner's correction also omits all 
hadron self-energy and vertex-renormalization diagrams. 
The assumption implicit in this approach is that such 
diagrams become part of what one calls the "electromag- 
netic form factor" of the struck hadron. 



VI. SIMULATION FRAMEWORK 

The models for (e, e'p) cross sections in PWIA, for ex- 
ternal bremsstrahlung, and for internal bremsstrahlung 
were implemented in the simulation code AEEXB |lq,E9[ . 
This code also includes facilities enabling a fairly com- 
plete simulation of experimental factors such as target ge- 
ometry, beam energy dispersion, ionization energy losses, 
and experimental acceptances. All these facilities were 
used in order to make the comparison as realistic as pos- 
sible. For practical reasons, certain classes of ionization 
energy losses were not included. Since our electron en- 
ergies are above the critical energy (for which radiative 
energy loss processes become more important than those 
due to atomic ionization), the ionization losses had a neg- 
ligible effect on our results. Ion-optical magnetic trans- 
port is also possible in AEEXB, using an interface to the 
standard ion-optics program TURTLE [ptpOl, but was 
neither necessary nor used for the current project. 

A complete cross-section simulation including the ra- 
diation tail consists of several distinct pieces which must 
be combined at the end to obtain the final result. The 
framework is sketched here; readers wishing to see a more 
detailed explanation should refer to |2]]] . The discussion 
below makes the simplifying assumption that the final 
state space for the residual (A— 1) nucleus consists of one 
discrete state plus a continuum; this condition is satis- 
fied for the reaction with which we compare, 3 He(e, e'p). 
Multiple discrete states would be straightforward to im- 
plement. The spectral function used for 3 He comes from 
the INFN/Rome group [||J2|]. 

A complete simulation consists of the following indi- 
vidual simulation runs: 

1. two-body breakup with external bremsstrahlung. 
This run handles computation of th e "un radiated" 
part of the cross section (see Sec. VD ). No in- 
ternal bremsstrahlung is computed; rather, the 
computed PWIA cross sections are reduced by 
the Schwinger correction to account for the frac- 
tion which will be redistributed into the internal- 
bremsstrahlung tail. Sampling is performed in 
(p e i,9 p ,4> p ); the constraint of a definite (A — 1) 
final-state mass provides the solution for p p < . Ex- 
ternal bremsstrahlung is allowed before the (e, e'p) 
vertex (modifying the beam energy) and after- 
wards (modifying the scattered-electron energy). 



The external-bremsstrahlung distribution is di- 
rectly sampled, obviating the need for a cutoff cor- 
rection. 

2. two-body breakup with internal and external 
bremsstrahlung. This run handles the part of the 
cross section which has been redistributed into 
the internal radiation tail. For each event, six 
variables are sampled (p e ',PpO- First, external 
bremsstrahlung is computed along the incident 
electron direction (possibly modifying the incident 
electron energy). Then solutions are found for 
the radiated-photon energies corresponding to in- 
ternal radiation along either the incident or scat- 
tered electron directions. The cross section is com- 
puted according to Eq. ([H]). Finally, external 
bremsstrahlung is computed along the scattered 
electron direction (possibly modifying the detected 
electron energy). Here only the virtual- photon part 
(1 — S v ) of the Schwinger correction is applied since 
the radiative tail corresponding to 5 r is what we 
are computing. 

3. continuum breakup with external bremsstrahlung. 
This piece is similar to case 1 above, except 
that events are sampled in six kinematic variables 
(Pe' , Pp' ) since there is a continuum of possible 
(A — 1) final states. 

4. continuum breakup with both internal and external 
bremsstrahlung. This simulation is similar to that 
of case 2 above, except that due to the complete 
kinematic freedom in the final state, the photon en- 
ergy is constrained only to be larger than AE (the 
experimental resolution). Therefore, seven vari- 
ables are sampled, the radiated-photon energy be- 
ing the seventh. 

5. detection volume simulation. This simulation is 
standard procedure for determining what fraction 
of the six-dimensional acceptance in (p e < , p p < ) can 
contribute to any given bin in a cross-section spec- 
trum. The results of this piece are used to properly 
normalize the simulated spectra when producing 
cross section results. No energy losses effects are 
included, since this part of the simulation only mea- 
sures the relative probability of detection of various 
kincmatical configurations (regardless of their ori- 
gin). 

The simulations, when properly weighted by sampling 
volumes and numbers of trials, are combined to form sim- 
ulated cross sections. The cross sections can be plotted 
as the same sort of spectra shown in experimental pa- 
pers, by sorting the simulated events into histograms in 
the same way an experimenter would sort data. 
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FIG. 5. Experimental and Theoretical Momentum Distri- 
butions for 3 He(e,e'p) 2 H. Experimental distributions are 
shown for three different electron beam energies; the figure 
of 0.84 in the text refers to the ratio of the 855 MeV experi- 
mental distribution to the theoretical distribution. 



VII. EXPERIMENTAL DATA 

The data with which we compare our simulations was 
acquired with the three-spectrometer detector setup at 
the MAMI accelerator facility in Mainz fHlf] . Two of 
these spectrometers were used to detect scattered elec- 
trons and knocked-out protons; the third served as a lu- 
minosity monitor. A cryogenic gas target provided the 
3 He target nuclei. The experiment measured cross sec- 
tions for the reaction 3 He(e, e'p) in a variety of kinematic 
settings. For more information on the experiment and its 
physics goals, the reader can consult Q; here we focus 
only on the essentials needed for the radiation-tail com- 
parison. The kinematical settings and experimental ac- 
ceptances for the data discussed here are given in Table 

B 

One question which must be addressed in this study is 
"how well can the PWIA model describe the reaction?" 
If the description is not favorable, further work is useless 
since our model computes the basic (e, e'p) cross section 
in PWIA. 

Figure H shows the measured experimental momentum 



distribution n B 



<T/{Kcr e 



kinematical factors in Eq. 



where K stands for the 
The experimental cross 



section a above has been radiatively corrected using the 
traditional technique, which has been shown to work well 
in this experiment [Q for excitation energies less than 
20 MeV. After radiative correction, the two-body peak 
could be cleanly resolved from the continuum and a cross 
section assignment is straightforward. 

The momentum distribution is compared to the 
theoretical two-body breakup spectral function J22 23 
S(S p ,p m ). S p is the single-proton separation energy and 
corresponds to 3 He — ► p + 2 H. Aside from an overall 
scaling factor of 0.84, the theoretical spectral function 
is in good agreement with the data. We interpret this 
agreement as an attenuation of the outgoing proton flux 
in the reaction, due to FSI, of constant magnitude 0.84; 
aside from this, effects outside the PWIA are not impor- 
tant. Ref. fm shows several other instances of how, apart 
from this overall reduction, PWIA calculations describe 
the data well. This good agreement can be attributed 
to our use of a light nucleus (reducing FSI effects) and 
the fact that our kinematics are directly tuned to the 
quasi-elastic point, where PWIA should work best. 



VIII. RESULTS 

We will discuss the results in two stages. First we 
will present results using the unmodified Boric-Drechsel 
tail computation (/ mp = 1), including our extension for 
discrete states of the (A — 1) residual nucleus. These 
results show a clear discrepancy in the tail region. We 
then present the derivation of our tail correction factor 
/ mp . Then we present results including this correction, 
which will be shown to resolve the discrepancy. 

Since the spectral function falls rapidly with p m , we 
were concerned about relying on the theoretical momen- 
tum distribution over the large p m acceptance of this 
experiment. Eventual discrepancies between our calcu- 
lation and the experiment might be due to inaccuracies 
in the hadronic structure of 3 He. In order to reduce this 
possibility, we carried out this study within a limited 
regime of p m by placing a cut on both the experimental 
data and on the simulation results. For all plots shown 
below, only missing momenta in the range 40 < p m < 50 
MeV/c are considered (for both the experimental data 
and the simulation). This particular region is near the 
top of the experimental acceptance in p m , where we make 
the best measurement of the two-body momentum distri- 
bution (on which the scaling factor is based) and where 
we had the greatest statistical accuracy in the experi- 
mental tail cross section. 



A. Results with Unmodified Tail Cross Section 

Fig. ra shows a comparison of the measured cross sec- 
tion, plotted as a function of the measured missing en- 
ergy E m , and the results of our simulation at the same 
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FIG. 6. Differential cross sections measured in the Mainz 
experiment (histogram) and computed by the simulation pro- 
gram (solid line). The experimental data have not been cor- 
rected for bremsstrahlung effects. The multi-photon tail cor- 
rection / mp has been set to unity. 



kinematics. The simulation result has been scaled by the 
factor 0.84 in accordance with the findings for the mo- 
mentum distribution. 

The agreement is generally excellent, the shape having 
been perfectly reproduced within the statistical accuracy 
of the simulation. The differences at the low-i*^ side 
of the peak in the spectrum are not really worrisome, 
since our simulation did not include all possible mecha- 
nisms of energy loss and its accompanying contribution 
to the experimental resolution. For example, while ex- 
ternal bremsstrahlung in the target-cylinder walls was 
accounted for, ionization energy losses in this material 
(82 /im foil of iron), for the incident and scattered elec- 
trons and ejected proton, were not. Thus sharp features 
in the cross section (such as the \ow-E m peak) will not 
be correctly reproduced by the simulation. At missing 
energies below about 10 MeV, where the radiative tail 
is still a small contribution, the integrals of the simula- 
tion and the data agree to within 1% (the integrals are 
not sensitive to the shape differences discussed above). 
This directly indicates that the empirical scaling factor 
is applicable to the continuum breakup as well, since our 
scaling factor of 0.84 was fixed by the two-body breakup 
results alone. For E m above 10 MeV, the simulation 
predicts a larger cross section than observed, with an es- 
sentially constant excess of about 20% (see Fig. 0). Since 
the shape reproduction is excellent, and the strength in 
the low-.E m region is well described, the comparison sug- 
gests a problem with the amplitude of the calculated tail, 
but not with its shape. 
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FIG. 7. Ratio of simulated (e, e'p) spectrum (normalized 
by 0.84 as discussed in the text) to the experimentally mea- 
sured cross section spectrum, plotted as a function of missing 
energy. A one-parameter fit to these data over the region 
10 < E m < 75 MeV yields a ratio of 1.18 (shown in the fig- 
ure). The bin-to-bin fluctuations in the ratio are due to the 
statistical uncertainties in both the cross section and the sim- 
ulation. The only significant deviation from the fit is below 
10 MeV, where the ratio decreases towards unity. 



B. Radiation Tail Correction for Multiple-Photon 
Processes 

The Schwinger correction, which has been applied to 
the "unradiated" strength dominating the region of low 
E m , includes the effects of multiple-photon emission and 
here the simulation and data agree. The tail region cross 
section has been derived to first order in real-photon 
emission, and here the simulation does not agree with 
the measured data. Multiple photon processes are there- 
fore clearly indicated as a likely source of the discrepancy 
in the tail strength. 

A rigorous derivation of a multiple-photon tail cross 
section is beyond the scope of this paper, but an in- 
tuitive derivation is easy to provide. In the limit that 
the variation in the PWIA vertex cross section is very 
slow, bremsstrahlung processes only redistribute strength 
with respect to the asymptotic kinematics. Thus if we 
add the "unradiated" part still residing in the peak to 
that residing in the tail, we should recover the original 
PWIA cross section cr pw i a . 5 r represents the fraction of 
strength radiated out of the peak, to first order. Thus if 
the Borie-Drechsel cross section is valid in first order, its 
integral will also yield a fraction S r of <7 pw i a . However, 
the fraction remaining in the peak has been adjusted to 
account for higher-order radiation; the cross section here 
is e _l5r CTp W i a . The sum of the two is (e _<5r + 5 r ) Cpwia- 
The factor in parentheses is differs from unity in second 
order. If we apply the multiplicative factor 
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FIG. 8. Differential cross sections measured in the Mainz 
experiment (histogram) and computed by the simulation pro- 
gram (solid line). The experimental data have not been cor- 
rected for bremsstrahlung effects. The full multi-photon tail 
correction / mp has been used. 
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FIG. 9. Ratio of simulated (e, e'p) spectrum (including the 
full tail correction and normalized by 0.84 as discussed in the 
text) to the experimentally measured cross section spectrum, 
plotted as a function of missing energy. A one-parameter fit to 
these data over the region 10 < E m < 75 MeV yields a ratio 
of 0.98 (shown in the figure). The bin-to-bin fluctuations in 
the ratio are due to the statistical uncertainties in both the 
cross section and the simulation. 



to the tail cross section, we recover o- pw i a for the 
sum of peak and tail cross sections in the presence of 
bremsstrahlung. 

Note that this discussion only concerns the real-photon 
part of the internal bremsstrahlung correction. The ex- 
ternal bremsstrahlung distribution described above is an 
energy-loss distribution which includes the higher-order 
contributions, thus they do not need to be considered 
here. 



here we consider an expanded range of p m . This check 
was made to ensure that our agreement had not been 
fine-tuned for only the small region of p m we had been 
considering. The corresponding ratio plot is very simi- 
lar to Fig. g, with the one-parameter fit yielding a ratio 
0.976. The same simulation produced both Figs. @ and 
|l0| ; the only difference was a change in the p m condition 
specified in the histogram-sorting program. 



C. Simulation including multi-photon tail correction 

The simulation was repeated including the multi- 
photon tail correction, but otherwise identical (includ- 
ing the scaling factor 0.84). The results are presented in 
Figs. H and [| 

The reproduction of the shape of the tail is still excel- 
lent, which is not surprising. For the chosen kinematics, 
S r has an average value of 0.46, with a la deviation of 
only 0.8% across the physical acceptance — / mp is essen- 
tially a multiplicative constant for the entire tail. How- 
ever, the simulation now reproduces the strength of the 
tail to the same level of accuracy as for the peak region. 
These excellent results provide unambiguous proof that 
multiple-photon processes are important in the radiation 
tail, and also that our proposed correction factor is valid 
at the few percent level (at least in the kinematical regime 
studied here). 

Finally, we show in Fig. 



10 a similar comparison of 



experimental and simulated cross-section spectra, except 



D. Decomposition of Cross Section 

It is instructive to separate this cross section calcula- 
tion into its components. This is a luxury that nature 
does not afford the experimenter. One such decomposi- 
tion is shown in Fig. [jl]. Recall that what one usually 
wants to measure is the cross section with the radia- 
tion tail removed. This corresponds to the dotted line 
in Fig. O, which is the simulation result for continuum 
breakup^He(e, e'p)np with bremsstrahlung turned off. 
The dashed curve shows the simulation for 3 He(e,e'p)d 
only, but including the full bremsstrahlung tail. The solid 
curve is the total simulation result as shown in Fig. |8j. 

It is immediately apparent from the figure that there is 
no hope of making a significant measurement for missing 
energies much above 15 MeV — statistical fluctuations 
associated with the tail subtraction procedure will render 
such a measurement insignificant. Most of the observed 
cross section for E m > 20 MeV is due to 3 He(e, e'p)d 
reactions residing in the radiation tail. A new experi- 
ment, made in a kinematical regime that does not result 
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FIG. 10. Comparison of experimental data to simulation 

(including tail correction). The plot shown is identical to that 

of Fig. H except for the current plot, events (both experimental 

and simulated) with missing momenta 30 < p m < 100 MeV / c 

are included. The bin size has been reduced to 0.25 MeV due 

to the much better statistical precision of these data. 
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FIG. 11. Decomposition of simulated (e, e'p) cross section. 
The solid curve gives the final result for the cross section 
including the entire spectral function and all radiative pro- 
cesses. The long-dashed curve only includes the two-body 
3 He(e, e'p)d part of the spectral function, along with both 
classes of bremsstrahlung. The dotted curve includes only 
the continuum 3 He(e, e'p)np part of the spectral function, and 
the bremsstrahlung effects are not included. When measur- 
ing in the continuum, the dotted curve is what one attempts 
to extract, and the long-dashed curve is physical background 
which must be removed by a radiative-correction procedure. 
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Missing Energy (MeV) 
FIG. 12. Experimental data for 3 He(e, e'p) before and after 
application of radiative corrections. In this figure, the exper- 
imental spectral function is displayed, which is related to the 
experimental cross sections by Eq. 



in the generation of such a strong radiation tail, will be 
required to make a statistically-significant measurement 
of the cross section in this region. 

We expect simulations such as that described here will 
become a standard tool in the planning of experiments 
at large missing energies, since one would clearly like to 
avoid performing an experiment in kinematical regimes in 
which the radiation tail is stronger than the cross section 
to be measured. 



IX. IMPLICATIONS FOR RADIATIVE 
CORRECTIONS OF (e,e'p) DATA 

We stated in the introduction that this work was begun 
in an effort to understand problems encountered when 
applying radiative corrections to the 3 He(e, e'p) data dis- 
cussed here. Fig. O (from Ref. O) illustrates the prob- 
lem. 

The figure compares the measured spectral function 
(from the same dataset which produced Fig. pi) with the 
corresponding spectral function after applying radiative 
corrections, i.e. removing the radiation tail. The tail cor- 
rection procedure is essentially identical to that of g] 
which was briefly described in Sec. IV. For E m > 25 



MeV, the corrected spectral function is negative, clearly 
indicating a deficiency. 

In this case, the defect can not be obviously traced to 
multiple-photon emission, since the tail computation is 
based on the distribution of the exponential form of the 
Schwinger correction in E m . Specifically, the number 
of expected experimental counts inside a certain missing 



p(0) 



energy bin E m < E m < Em + AE m is given by 



(o) 



N cxp (E£\AE m ) = e- 



S r (AE„ 



] N (EM) 



Here TVn is the number of counts which would have 
been measured in the absence of bremsstrahlung, and 
E m is the missing energy at which the counts would 
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have appeared in the absence of bremsstrahlung (or other 
processes which modify the asymptotic particle energies 
such as ionization energy loss). As AE m becomes larger, 
the bin includes more of the radiation tail (the unra- 
diated strength is included by definition) so that iV cxp 
approaches Nq. The distribution of the radiation tail in 
E m is thus 



dN c 



dE n 



N, 



d(e 



-M 



d(AE m ) 



The factor e~ Sr again explicitly includes the multiple- 
photon processes. 

There are, however, several possible other reasons for 
the failure of the radiation-correction procedure. We dis- 
cuss them in detail in the following sections. 



A. Incomplete Kinematic Reconstruction 

The radiative-correction process is usually applied in 
two kinematic dimensions (see Ref. [nl for a thorough 
discussion). A common choice for the two dimensions 
is (E m ,p m ). First, a two dimensional cross-section his- 
togram is created with the independent variables being 
E m and p m . Fig. would be appropriate if the z-axis cor- 
responded to cross section. As discussed in Sec. jM the 
correction would begin at the left-hand edge of Fig7|| For 
each bin in p m at this bin in E m , correction factors would 
be applied and tails would be generated and subtracted 
for the bins to the right. In such a procedure, there is 
no information about any of the other kinematic param- 
eters; a swath in seven-dimensional kinematic space has 
been reduced to a two-dimensional pixel. The most com- 
mon remedy for this lack of information is to treat the 
entire bin as if the rest of the parameters were fixed at 
their central values. However, Fig. shows a substan- 
tial dependence in the tail trajectories on the relative 
angle between q and p m . While both of these may vary 
across the detector acceptances, both are held fixed in 
the correction procedure. |q| is also held fixed, while we 
know that it varies substantially across the acceptances 
and causes large changes in the cross section through 
a ep . For example, at the kinematics corresponding to 
Fig. ||, KQ 2 )/KQ 2 ) (the standard deviation of Q 2 di- 
vided by the mean value) is about 7.6% for the pixel 
5.5 < E m < 5.9 MeV, 40 < p m < 50 MeV/c. This leads 
to a 15% RMS variation of the cross section due to o~ ep 
alone. The variation in the computed 3 He(e,e'p) cross 
section is about 20%. 

Procedures have been developed to perform the correc- 
tion in more dimensions (e.g. four were used in Ref. Jlj]) 
but such schemes are only feasible for experiments with 
good statistical precision, as each additional dimen- 
sion tends to reduce the statistical precision per bin by 
roughly an order of magnitude. 



B. Simplifying Assumptions About the Radiation 
Tail 



The standard correction procedure uses a derivative of 
the Schwinger correction factor vs. E m to generate the 
tail distribution. However, as explained in Sees. E\] and 
M, there are two directions this tail can take in (E m ,p m ) 
space, corresponding to the two terms of the peaking 
approximation. The Schwinger correction gives no guid- 
ance as to how much strength resides in each tail. The 
standard practice is to assume [| 12 1: 



• the incoming and outgoing electrons contribute in- 
dependently to the tail, so one may factor e~ s = 

C e (AE m )C e ,(AE m ); 

• the two contributions are equal, so C e = C e > = 

The Borie-Drechsel formula (Eq. dl4)) for the radiation 
tail clearly does not have these properties. Firstly, the 
two tails add instead of multiply. Secondly, they are not 
equal. Even for vanishingly small photon energies (k e 
and k e >), the two terms differ by the factors \n(2E e /m e ) 
vs. In (2E e > /m e ) (which have values 8.12 and 7.81 respec- 
tively in the kinematics studied here). As the radiated- 
photon energy increases, the difference between the two 
terms also increases; for a photon energy of 100 MeV, 
the incident-electron contribution is 6% larger than that 
of the scattered electron in the present kinematics. The 
difference between the tail magnitudes is mainly driven 
by the ratio tv/E e ; when it is large, the tail strengths 
differ more. For a specific experiment planned at JLab 
with u) = 834 MeV and E e = 1245 MeV, the two tails 
differ by about 16% in strength. 



C. Comparison With Direct Tail Simulation 

The radiation tail calculation presented here suffers 
from none of the above deficiencies. 

• the tail is generated event-by-event, so for each tail 
evaluation, the complete kinematic information is 
available 

• the distribution of the tail strength in this kine- 
matic space is based on a first-order QED calcula- 
tion, not on plausible assumptions. 

The main deficiency of our computation is the nature of 
the multi-photon correction factor. It is based on argu- 
ments of probability conservation rather than on a rigor- 
ous QED calculation. This argument is however of the 
same type which leads to the exponentiation of 5 r in the 
standard approach. A critical review of including higher- 
order terms via exponentiation (including a summary of 
relevant literature) can be found in [|16| . 
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D. An Improved Radiative Correction Procedure 

The findings reported above indicate that the stan- 
dard radiative correction procedure is not likely to work 
for cases in which the detector acceptances are relatively 
large (producing large variations in Q 2 for individual pix- 
els in cross-section histograms) or in which a correction 
is being made over a large range in E m (so that the dif- 
ferences between the two tails becomes important). Our 
findings suggest an improved method for radiatively "cor- 
recting" experimental data. 

The procedure would begin with a model spectral func- 
tion and an accurate model of the experimental appa- 
ratus, such as has been described here. A simulation 
code similar to ours should be used to generate a "ra- 
diated" cross section spectrum. A comparison between 
the experimental and simulated histograms will indicate 
regions of discrepancy, and the discrepancy function can 
be used to modify the model spectral function. The pro- 
cedure is repeated until it converges, at which point the 
model spectral function corresponds to the unradiated 
result. This procedure is independent of the PWIA if we 
replace the theoretical spectral function described above 
with the "distorted" spectral function |26|]. We learned 
during the final stages of preparing this article that such 
an iterative procedure has been developed and success- 
fully applied for an experiment in Hall C at Jefferson Lab 



X. FURTHER WORK 

It is desirable to have a more rigorous theory provide 
the multi-photon tail cross section. The beginnings of 
such an approach can be found in Q . We encourage this 
group to complete and publish these results, especially 
since they have made some detailed evaluations of their 
approach in the Jefferson Lab energy domain. 

It would also be interesting to reanalyze some of the 
older high-ii'm (e, e'p) data using the improved technique 
described here, since these data have been a source of 
controversy, given the sometimes puzzling behavior of the 
cross section at high missing energies. The current study 
indicates that some of this cross section might well be 
misidentified bremsstrahlung strength. 



XI. CONCLUSIONS 

We have presented a framework for computing (e, e'p) 
cross sections which includes the radiation tail to first 
order. The computed cross sections have been compared 
to experimental data in such a way that effects such as 
acceptance averaging are correctly accounted for, allow- 
ing a direct evaluation of the radiation tail cross section 
calculation. 



The computed tail reproduces both the shape and 
magnitude of the experimental spectrum perfectly within 
experimental errors. It was necessary to derive a 
correction, applied to the radiation tail, for higher- 
order bremsstrahlung effects before this agreement 
could be obtained; the original tail calculation treated 
bremsstrahlung only to first order. 

A standard radiative correction procedure has also 
been applied to these data. Such a procedure is designed 
to move the radiation tail strength back into the origi- 
nating kinematic bins. The straightforward application 
of this procedure (that is, without tweaking parameters 
to improve the agreement) results in physically unrea- 
sonable "deradiated" cross sections in the tail-dominated 
part of the spectrum. There arc several reasons to expect 
such a failure. An obvious one is that this procedure 
collapses a complicated kinematical hypersurface (along 
which the cross section varies substantially) to a single 
point in (E m ^p m ). More subtle are the disturbing differ- 
ences between the properties of the correction-procedure 
tails and those of a tail cross section rigorously computed 
in QED. The observed flaws in the correction procedure 
are not likely to affect earlier data taken at low missing 
energies, e.g., at NIKHEF, Bates, and Mainz. They may 
affect earlier high-E^ data, and will likely be fatal for 
several of the (e, e'p) experiments planning to measure 
at large--E m at Jefferson Lab. Our results indicate how 
radiative corrections should be applied so as to avoid such 
problems. 

The current project has yielded quantitative illustra- 
tions of the failure of the standard radiative-correction 
procedure for (e, e'p) experiments. We have also shown 
that a simulation, coupled with an accurate model for 
the radiative-tail cross section, can radiate the theory 
(instead of deradiating the data) and achieve excellent 
agreement with experiment in a situation where the 
correction procedure fails. The simulation technique 
described here provides a basis for iterative radiation- 
correction procedures for future (e, e'p) experiments. 



APPENDIX A: OTHER LIMITATIONS OF THE 
PWIA 



In Sec. Ill we were careful to distinguish the spectral- 



function quantities e and p from the experimentally de- 
termined values E m and p m . Even in the absence of 
bremsstrahlung this is necessary since the PWIA never 
holds completely, and is sometimes grossly violated. For 
such cases, e ^ E m and p ^ — p m . Final State Interac- 
tions (FSI) between the ejected proton and the residual 
nucleus provide an illustrative example of how the cor- 
respondence is broken. At the photon-proton vertex of 
Fig. EL the amplitude for the interaction will depend on 
the particular values of e and p. A subsequent interaction 
between the ejected proton and the residual nucleus can 
change both the momentum and excitation energy of the 
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residual system, thus leading (through Eq. (g)) to values 
for E m and p m different than e and p. While this point 
is not particularly relevant for the present work, we men- 
tion it here for completeness and because it is apparently 
often overlooked. 
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